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Q ! ABSTRACT 

C/3 , Concatenating data from the millimetre regime to the infrared, we have performed 

■ sp ectral energy dist ribution modelling for 227 of the 405 millimetre continuum sources 

' ' of lHill et all (I200I which are thought to contain young massive stars in the earliest 



stages of their formation. Three main parameters are extracted from the fits: tem- 
^ ■ perature, mass and luminosity. The method employed was Bayesian inference, which 

00 ' allows a statistically probable range of suitable values for each parameter to be drawn 

10 ■ for each individual protostellar candidate. This is the first application of this method 

to massive star formation. 

CO I The cumulative distribution plots of the SED modelled parameters in this work 

. indicate that collectively, the sources without methanol maser and/or radio continuum 

associations (MM-only cores) display similar characteristics to those of high mass star 
00 , formation regions. Attributing significance to the marginal distinctions between the 

■ MM-only cores and the high-mass star formation sample we draw hypotheses regarding 

the nature of the MM-only cores, including the possibility that the population itself is 
comprised of different types of source, and discuss their role in the formation scenarios 
' of massive star formation. In addition, we discuss the usefulness and limitations of SED 

^ I modelling and its application to the field. From this work, it is clear that within the 

valid parameter ranges, SEDs utilising current far-infrared data can not be used to 
determine the evolution of massive protostars or massive young stellar objects. 

Key words: stars: formation ~ stars: fundamental parameters - stars: early-type - 
submillimetre - masers - (ISM:) Hii regions. 



1 INTRODUCTION 

The study of massive stars and their natal molecular envi- 
ronments has seen a surge in interest and study in recent 
years. Despite this, precise scenarios for the formation of 
massive stars, especially of the earliest stages of their evolu- 
tion, are still not forthcoming. This may be attributed to a 
number of factors. The very nature of massive star formation 
(i.e. rare, rapid, clustered and distant) impedes study of the 
earliest stages of their evolution, which are not easily distin- 
guished. Massive stars also form in turbulent and evolving 
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environments, e.g. jets, disks and outflows, which contribute 
to hinder extraction and interpretation of information from 
these regions. Additionally, instrumental limitations, specif- 
ically the resolution of current instruments, are not suflicient 
to probe the inner most workings of the cocoons in which 
massive stars are forming. 

Previous studies have identified associations be- 
twce n young massive stars and methanol masers 
(Pes talozzi et all [2005il. Uhra Compact Hn (UC Hn) 
regions (iThompson et a l. 2006), IRAS colour selected 
sources (|Wood fc ChurchwelL .19891) as well as MSX colour 



20021 ). Methanol masers 



selected sources (jLumsden et al. 
and UC Hii regions, in particular, are thought to fea- 
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ture prominentl y in the ear l iest s t ages of massive star 
formation (c.f. iBatrla et al. 1987: ICaswell et al. 1 1995 



n- 



995| : 
003; 



Minier et al. 2001; Bcuthcr ct al. 2002; Faiindez et al. 2004 
Williams o t al. 2004; Thompson ct al. 2006, and references 



within) . 

Early work (e.g. IWalsh et al.lll998l l suggested that the 
methanol maser was the earliest indicator of massive star 
formation prior to the onset of Hii regions that are sign- 
posted by radio continuum emission. More recent work has 
focused on finding a precursor stage to that signposted by 
radio continuum emission, which would mark the very earli- 
est stages of massive star formation. The hot molecular core 
(HMC) and infrared dark cl ouds are such objects proposed 
to satisfy this criter i on (cf. Olmi et al. 19961 ; lOsorio et al.l 



1 19991 : iHill et aLllioOSI : iRathborne et alJl200'^') The me t hanol 
masers likely form during these stages " ifHiU et all l2005l : 
iLongm orc ct al. 2006, 2007) 

[kill et al.l (|2005 . hereafter Paper I) undertook a SIMBA 
millimetre continuum emission survey toward regions dis- 
playing evidence of massive star formation, in search of cold 
cores that would mark the earliest stages of their evolution. 
This survey revealed each of the methanol maser and radio 
continuum sources targeted to be associated with millimetre 
continuum emission. Interestingly, this survey also revealed 
evidence of star formation clearly offset from, and devoid 
of, both the methanol maser and radio continuum sources 
targeted. These sources were dubbed 'MM-only' corefl. 

Preliminary analysis showed these MM-only sources 
to be smaller and less massive than cores harbouring a 
methanol maser and/or an UC Hii region. This conclusion 
was drawn however assuming a constant temperature of 20 K 
across the sources in the sample. At least 45 per cent of 
these 'MM-only' sources are also without mid-infrared MSX 
emission. It was consequently proposed that the MM-only 
core is a possible precursor to the methanol maser stage 
of massive star formation, and thus traces an even earlier 
stage - perhaps even the earliest stage in the formation of 
massive stars. F ollow-up submil limetre observations of these 
MM-only cores l|HiU et al.ll2006l . hereafter Paper II) revealed 
each of them to be associated with submillimetre continuum 
emission confirming their association with cold, deeply em- 
bedded objects. 

In order to ascertain the nature of the MM-only cores 
and any role that they play in the formation and evolu- 
tion of massive stars, it is necessary to determine their am- 
bient physical conditions, such as temperature, luminosity 
and mass. Only in light of this information is it possible to 
characterize the MM-only core and address hypotheses re- 
garding their formation and/or whether they are indicative 
of the earliest stages of massive star formation. Addition- 
ally, in order to put the MM-only core into context within 
an evolutionary sequence for massive star formation, a wide 
cross-section of sources suspected of being at various evolu- 
tionary stages is required. 



^ In this paper, the term 'cores' refers to molecular cloud frag- 
ments that were detected through millimetre dust continuum 
emission. These 'cores' have sizes and masses that span a large 
range which qualifies them to form many young stellar objects or 
even protostellar clusters. Although they are more widely labelled 
as 'clumps' in the literature, for consistency with our earlier work 
we preferentially use the term 'cores' in this work. 



In this paper we combine the (sub) millimetre data from 
our earlier wo rk (Papers I and II), w i th existing submil- 
Umetre data (iPierce-Price et al.l l2000l : IWalsh et all l2003l : 
[Thompson et ah 20061 ). together with archival MSX data 
where applicable, and IRAS data (more often than not as 
upper limits) to draw spectral energy distribution (SED) di- 
agrams for a large sample of sources detected with SIMBA. 
The sample itself is a cross-section of sources suspected of 
representing different evolutionary stages of massive star for- 
mation (see ii l2.ip . 



2 MULTIWAVELENGTH COMPILATION OF 
THE SIMBA SOURCES 

In this section, we outline the source selection criteria, ex- 
plore each of the wavebands used for the SED fitting, and 
explore infrared associations for the sample. 



2.1 The sample origin: (sub)millimetre 
observations 

Our 1.2 mm SIMBA survey (Paper I) revealed a total of 405 
millimetre continuum sources, a large number of which are 
MM-only cores as introduced in section [1] 

The SIMBA sample is comprised of four distinct classes 
of source, distinguished by the presence, or lack thereof, 
of methanol maser and radio continuum tracers. Class M 
sources are millimetre sources with methanol maser sites 
but are devoid of radio continuum emission. Class R sources 
are millimetre sources with radio continuum emission but 
without methanol maser emission. Class MR sources are 
millimetre sources with both methanol maser and radio con- 
tinuum emission. The fourth class of source is the MM-only 
sample which is comprised of sources with millimetre con- 
tinuum emission, but without methanol maser sites or UC 
Hii regions. 

The observation and data reduction method of each of 
the SIMBA and SCUBA surveys, are as described in Papers 
I and II, respectively and the reader is referred to these for 
more information. 

The millimetre (1.2 mm) and submillimetre (450 & 
850/im) fluxes used in this work were extracted directly from 
the (sub)millimetre continuum maps using the respective 
data reduction and analysis packages (see Papers I and II). 
This procedure involved distinguishing the source from the 
background and subtracting the latter from the former using 
apertures defined to 10 per cent contour level (of the peak 
flux) for each of the SIMBA and SCUBA sources, respec- 
tively. 

The definition of the source size (contour level of 10 
per cent) is critical for SED analysis as it infiuences the 
amount of fiux input into the SED as well as the resultant 
parameters from the model. Comparison of the integrated 
fiux determined for a sample of sources to a contour level of 
5, 10 and 20 per cent of the peak flux of the source, reveals 
an integrated flux difference of less than 10 per cent for a 5 
per cent contour and less than 15 per cent for a 20 per cent 
contour, when compared with our assumed size of a 10 per 
cent contour. As discussed in Section 13.11 we assume a 20 
per cent flux error for SED analysis which accounts for any 
ambiguity in source size. 
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2.2 Infrared data 

2.2.1 Infrared Astronomical Satellite (IRAS) 

The Infrared Astronomical Satellite (IRAS) was a joint sci- 
entific project between the United Kingdom, the Nether- 
lands and the United States. IRAS' mission was to per- 
form a sensitive and unbiased all-sky survey centred at four 
wavebands in the infrared regime: 12, 25, 60, 100 fim. IRAS 
was launched in January 1983 and ended its mission ten 
months later in November, after surveying 96 per cent of 
the skjjfj. The angular resolution of IRAS varies between 
about 30 arcsec at 12fim to about 2 arcmin at 100/im. For 
m ore detail about the des ign and performance of IRAS refer 
to lBeichman et all (119881 ). 

The IRAS Sky Survey Atlas (ISSA) is a publicly avail- 
able set of FITS images of the infrared sky, while the IRAS 
point source catalogue (PSC) provides flux estimates of 
sources detected by IRAS for each of the wavebands of ob- 
servation. 

The IRAS infrared fluxes at 60 and 100 /im have been 
extracted from the IRAS PSC for those sources with an 
IRAS association. They are also used for the purposes of 
upper limit constraints to spectral energy distribution anal- 
ysis in the absence of a direct association. Note that we do 
not use the other two IRAS wavebands at 12 and 25 ^m 
as they overlap that of the MSX satellite, which achieved 
better angular resolution (see following section [2.2.2p . 

2.2.2 Midcourse Space Experiment (MSX) 

The Midcourse Space Experiment (MSX) was a multi- 
discipline experiment sponsored by the Ballistic Missile 
Defence Organisation. Launched April 24, 1996 this in- 
frared satellite operated at a temperature of 11 to 12 K, 
and spanned the infrared regime from 4.2- 26 /^m. The four 
main wavebands of MSX are centred at 8.3, 12.1, 14.7 and 
21.3 fim. 

The MSX Galactic Plane Survey mapped the Galactic 
Plane for |6| < 5°, and surveyed the part of the sky missed by 
IRAS in the "Survey of Areas Missed by IRAS' , as well as 
other surveys. An overview of the astronomic al experiments 
conducted with MSX is given bv lPricel(|l995l ). while a com- 
plete de scription of t he ex periments and data processing is 
given in lPrice et all l|200ll) . 

The infrared instrument on MSX SPIRIT III, had a 
spatial resolution of 18.3 arcsec, and a sensitivity of O.lJy 
at 8.3fim. 

The MSX images were examined for sources appearing 
in the SIMBA source list. For those sources with a direct 
MSX association, the mid-infrared flux density was deter- 
mined from the calibrated images. The maps were converted 
from B1950 Galactic coordinates to J2000 equatorial coor- 
dinates. The maps were then converted from Wm^sr"^ to 
Jy, allowing for 6 arcsec square pixels, and an additional 
factor of 1.133 to convert from square pixels into a gaussian 
area. The final conversion factors for each of the wavebands 
were 6.84xl0^ 2.74x10*, 3.08x10* and 2.37x10* Jy per 
Wm'^sr"^ at 8.3, 12.1, 14.7 and 21.3 /xm, respectively. 

The flux of the source was then measured using the 

^ See |http: / /lrsa.ipac.caItech.edu/IRASdocs/iras.html| 



karma/kvi^ package by applying an aperture around the 
sources and at various points in the image considered to be 
the background. Contour levels of 10 per cent of the peak 
source flux were overlaid and used to define the 'source' aper- 
ture size. The resultant flux of the source was then deter- 
mined by measuring the flux inside each of the source and 
background apertures and subtracting the latter from the 
former. In most instances, the fluxes which were determined 
from the MSX images were consistent with (to within 10 per 
cent of) the flux reported by the MSX catalogue. However, 
as there are known problems with the fluxes reported in the 
MSX catalogue, we cautiously opted to manually determine 
fluxes from the MSX images for all sources. 

For those sources with an MSX association, the infrared 
MSX emission is used to draw the SED as described in sec- 
tion [3] 

2.3 Correlating the SIMBA data with other data 

The aim of this paper is to draw spectral energy distribution 
(SED) diagrams for individual sources in the SIMBA source 
list of Paper I. Submillimetre associations with these data 
were explored in Paper II. For those SIMBA cores which 
are resolved by SCUBA into multiple components, the cor- 
responding submillimetre flux for the SIMBA sources has 
been used to draw the SED, with individual submillimetre 
components summed together. 

Extracting source specific information from all sky sur- 
veys can be difficult due to confusion in the images. The in- 
frared images from the IRAS and MSX satellites are exam- 
ined here with respect to each of the sources in the SIMBA 
images. 

As a consequence of the poor spatial resolution of IRAS 
compared to that of SIMBA (24 arcsec), it was often not 
possible to conclude an IRAS association with an individ- 
ual SIMBA source. Generally the entire SIMBA map (240 
X 480 arcsec'^) falls within a single IRAS source, that is, the 
SIMBA instrument resolved the corresponding IRAS source 
into multiple millimetre components. Very few sources in 
the sample had a direct correlation with the IRAS peak of 
emission and hence an IRAS source. An IRAS fiux was gen- 
erally only used for isolated SIMBA sources. For all other 
sources, the IRAS flux corresponding to the nearest coin- 
cident methanol maser and/or radio continuum source was 
used as an upper limit in the SEDs. In both of these in- 
stances, the IRAS flux was taken directly from the PSC. In 
the few cases where a SIMBA source was completely devoid 
of IRAS emission, falling in diffuse background emission in- 
stead, an upper limit was obtained from the IRAS images 
according to the procedure described for the MSX data (sec- 
tion [2]2]2]), and is included on the SED plots. 

Considering the resolution of IRAS the usefulness of 
these data for characterisation and constraint of the cold 
component of our SED fits is questionable. It may be argued 
instead that Spitzer data would be a better choice, especially 
the 70 /im data which easily supersedes the resolution of the 
60 /.im IRAS data. We have however, elected not to use the 
Spitzer data for our SED fits primarily as these data are an, 
as yet, unpublished data set. While the data are available for 

^ http:/ /www. atnf.csiro.au/computing/software/karma/ 
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download, it contains many artefacts and saturated values, 
which we are unable to quantify. 

As mentioned in section 12.2.21 the MSX satellite pro- 
duced higher resolution images than that of IRAS. However 
as a result of an excess of mid-infrared emission in the fields 
examined, it is often not possible to distinguish individual 
associations due to confusion. For instance, there can be ex- 
tended PAH emission in the MSX 8, 12 and 14 /im maps, 
which makes it difficult to extract the fiux associated with 
dust emission from a core. Consequently, although a SIMBA 
source may not have a direct association with a MSX source, 
it does not eliminate that same source from having asso- 
ciated mid-infrared emission. Often, a SIMBA source falls 
within diffuse mid-infrared MSX emission, yet they are not 
directly associated with a MSX peak of emission (i.e. a MSX 
source). For any SIMBA source where a MSX association 
was ambiguous, MSX data were not used. MSX emission 
associated with multiple SIMBA sources was not used un- 
less it was obvious which SIMBA source is dominated by the 
MSX fiux. 

In the case of G 49.49-0.37, the SCUBA images reveal 
it to be quite complicated at both 450 and 850/im, resolving 
the SIMBA sources into multiple components. The SIMBA 
sources also are only partially sampled by our SCUBA data. 
Due to this inadequate sampling, as well as confusion for 
associations with MSX and IRAS, SED analysis was not 
performed for this entire region (20 sources). 



3 SPECTRAL ENERGY DISTRIBUTION 
ANALYSIS 

The concatenation of the (sub) millimetre and infrared data 
described in Section [2] enables the spectral energy distribu- 
tion of each of the SIMBA sources to be drawn. 



3.1 Modelling procedure 

We have modelled our sample according to a simple two- 
component model denoting a central warm core surrounded 
by a colder dust envelope (see iMinier ilaDllooi). The 'hot' 
component of this model is assumed to radiate as a black- 
body sphere, whilst the 'cold' component accounts for opti- 
cally thin emission from the dust. The emerging spectrum 
is then defined by: 

= [7rB,(rhot) RLt + B^Tcoid) M.oid k{u)] /(f (1) 

where is the fiux density of the source, Rhot is the radius 
of the hot component of the source, B^, is the Planck function 
for a temperature of Thot and Tcoid, M^oid is the mass of the 
cold component, d the distance to the source and «;(;/) is the 
mass absorption coefficient. In this instance, «:(;/) is assumed 
to vary as a power-law with k(A) — kq (A/Ao)~'^cm^ per 
gram of dust at Ap = 1.2 mm , as per the opacity mo dels of 
lOssenkopf fc Henninid (|l994l ) (c.f. lMinier et al.ll2005h where 
Ko = 1.0 cm'^ g~^. We assume a dust to gas ratio of 100. 

Equation [1] solves for four parameters: iihot, Thot, Tcom 
and Afcoid. Due to the ambiguities with SED fitting, it is 
not possible to constrain each of these four parameters in- 
dependently. The robust estimation for the range of validity 
of the parameters instead requires the potential correlations 



Table 1. Range of values explored for each of the parameters of 
interest from the fitting procedure using log-space sampUng. 



Parameter 




range 


no. of values sampled 


^cold 


2.73 


... 100 K 


100 


Mcold 


1 . . 


. 10" Mq 


100 


That 


100 , 


. . . 1500K 


25 


Rhot 


10^ . 


. . 10*^ Rq 


25 



between each of the parameters to be taken into considera- 
tion. With this in mind, we systematically explored a grid 
of models by varying all four free parameters in order to 
ascertain a solution. The range of values explored with the 
fitting procedure are summarized in Table [T] Spectral en- 
ergy distributions were calculated for each combination of 
these parameters, resulting in a total of 6.25 million syn- 
thetic SEDs. 

Comparisons between the models and the observations 
were drawn according to a reduced value. For each of the 
individual sources, we computed a table of reduced sam- 
pling of the whole parameter space. That is, the observed 
SEDs were compared to all synthetic SEDs according to the 
calculated reduced values. These tables of reduced x^ val- 
ues are used in the following section to estimate the range of 
validity for each of the different parameters via the Bayesian 
inference method. 

For the purpose of the x'^ calculations, each observa- 
tional fiux was assumed to have a fiux error of 20 per 
cent, except for those 450 /xm submillimetre and IRAS data 
where a 40 per cent error was used. The error estimates 
for the (sub) millimetre data are consistent with Paper II, 
whil st for the infrared d ata these estimates are consistent 
with iMinier et all (|2005h and deemed reasonable for fitting 
purposes. 

For those sources devoid of coincident mid-infrared 
MSX emission, we have fit the (sub)millimetre data for the 
col d component of Equ ation[T]only, similar to what was done 
bv IMinier et al.l ()2005f ). In the absence of an IRAS associ- 
ation we have used the IRAS flux as an upper limit in the 
fitting procedure. 

The free-free emission contribution to the millimetre 
fluxes of the sources in our sample is expected to be mi- 
nor. Comparison of cm-band radio continuum fluxes with 
our 1.2 mm continuum fluxes indicates that the free- free con- 
tamination is typically a few percent, even for bright sources 
such as G 5.89 -0.39, which has a free-free contamination to 
the 1.2 mm fiux of not more than 5 per cent. We therefore 
do not expect free-free emission to significantly infiuence our 
fits and do not consider it in the fitting procedure. 

The luminosity of each source was determined through 
integration of Equation[l] We opt to limit the range for lumi- 
nosity integration between 1.2 mm and 8 /im, corresponding 
to the region encompassed by our data. At longer wave- 
lengths, there are no data to constrain the SED curve, the 
shape of which is then reliant upon the model alone. For con- 
sistency, we thus also limit the shorter wavelengths to the 
actual data. We stress that the resultant luminosity is rep- 
resentative of the luminosity between the range integrated 
(1.2 mm and 8 fim) only and is not a direct representation of 
the bolometric luminosity of the source (i.e. it is a lower limit 
to the luminosity). In order to ascertain the effect that lim- 
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iting the integration range had on the resultant luminosity, 
we also integrated the luminosity over the wavelength range 
0.1 /im to 3 mm. Although the luminosity for this range was 
roughly only 5 per cent greater on average than the range 
encompassing our data (1.2 mm to 8 /xm), individual sources 
could vary as much as 40 per cent. Despite the differences 
in the luminosity determined for each range, varying the 
integration range had negligible effect on the shape of the 
cumulative distributions in Fig. [21 



3.2 Validity range of parameters 

To determine the range of validity for each of the four free 
parameters, or two free parameters for the single component 
analysis, we used a B a yesian inference method (IPress et al.l 



analysis, we used a ij a yesian mterence metnod l|-rress et al 
11993 : iLav et al.|[l997l : [Pinte et al.ll2007l : iPinte et al.ll2008l ) 



This technique allows us to estimate the probability of oc- 
currence of each parameter value. The relative probability 
of a single point of the parameter space {i.e. one model) 
is proportional to exp(— x^/2), where refers to the re- 
duced of the corresponding model. All probabilities are 
normalized at the end of the procedure so that the sum of 
the probabilities of all models over the entire grid is equal 
to 1. 

The Bayesian inference method of SED modelling does 
not only give the best value for each parameter, it also pro- 
duces a range of suitable values for each parameter, with 
each value having an associated probability of occurrence. 

The Bayesian method relies on a priori probabilities for 
the parameters. For the purposes of our analysis, we assume 
that we do not have any preliminary available information, 
choosing instead a uniform a priori probability, which cor- 
responds to a logarithmic sampling of the parameters. 

Figure [1] presents the best fit SED (i.e. the model with 
the smallest value) and the relative figures of merit (prob- 
ability distribution diagrams) estimated from the Bayesian 
inference method for the temperature and mass of the cold 
component of the fit for three individual sources. These re- 
sults were obtained from marginalization {i.e. summing) of 
the probabilities of all models, where one parameter is fixed 
successively to its different values. The resulting histograms 
indicate the probability that a parameter takes a certain 
value, given the data and assumptions of our modelling. The 
width of the probability curves is a strong indicator of how 
well the data are constrained (see Section I4.1|l . The radius 
and temperature of the hot component are considered as 
"nuisance parameters" here. That is, they are parameters 
that have an influence on the data but do not have a direct 
physical interpretation and are thus not of prime interest to 
us in this study. 

For each parameter 6 with a density of probability p{6), 
this range of validity is defined as the interval [^i , 62] where: 



p(ei)=p(6»2) and 



p{e)de = 7 



(2) 



with 7 = 0.68. The interval [9i, 82] is a 68 per cent confidence 
interval and corresponds to the 1 a interval for a Gaussian 
distribution of probability. Table[2]gives the range of validity 
of the parameters for the cold component of the SED fit, 
with min and max representing the lower and upper values 
of this range, respectively. 



Table 3. Table indicating how many of each class of source sat- 
isfy the different types of SEDs applied to the sample. Column 
1 indicates whether a two-component or single-component SED 
has been applied. Column 2 indicates the class of the millimetre 
source. Column 3 indicates the number of sources with an IRAS 
association, whilst column 4 indicates the number of sources that 
have been fit with the IRAS upper limits. 







Source 


IRAS 


Upper 






Class 


Fit 


Limit Fit 


(1) 




(2) 


(3) 


(4) 


TWO- 


(MM) 


MM-only 


3 


11 


COMPONENT 


(M) 


maser 


16 


15 


(88) 


(MR) 


maser+radio 


4 


16 




(R) 


radio 


7 


16 


SINGLE (92) 


(MM) 


MM-only 


5 


56 




(M) 


maser 


1 


18 




(MR) 


maser+radio 


1 


8 




(R) 


radio 





3 



4 ANALYSIS 

In this section we discuss the method of our fitting approach 
and analyse each of the three parameters derived from the 
SED analysis: temperature, mass and luminosity with re- 
spect to the different classes of source in the sample. 



4.1 Method 

Spectral energy distribution diagrams have been drawn for 
227 of the 405 sources appearing in the SIMBA list of Pa- 
per I (56 per cent). Of these, 135 are two-component SEDs 
modelling a hot and a cold component for a source, whilst 
the remaining 92 are single component fits which model the 
cold component of a source only. 

For 47 two-component SEDs, the cold component of the 
curve is ill-constrained as a consequence of the poor data 
sampling: only the 1.2 mm SIMBA flux with IRAS upper 
limits. In this instance the number of free parameters (TcoM 
& Mcoid) for the cold component exceeds the number of data 
points available for analysis and it is not possible to define a 
value nor constrain the resultant parameters from these 
fits. These 47 'bad' SEDs are consequently excluded from the 
following analysis and discussion, bringing the total number 
of SEDs analysed to 180 (44 per cent of the sample) and the 
number of two-component fits to 88. 

The three main parameters that we extract from the 
SED analysis are the dust temperature, dust mass and the 
luminosity. The 180 millimetre continuum sources comprise 
four distinct classes of source, as discussed in section 12.11 
Table [3] indicates how many of each class of source satisfy 
each of the single and two-component SEDs, as well as how 
many sources have a direct IRAS association. Those sources 
that were fitted with IRAS upper limits are identified by an 
" in Table H 

For each of the sources with a near-far distance ambi- 
guity, the near-distance value of each parameter is assumed 
in the following analysis (and for comparison with Paper I). 
As a check, in Paper I, we examined the results for sources 
with no distance ambiguity with the results from the sample, 
assuming the near distance for 197 sources with a distance 
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Table 2. Parameters resulting from spectral energy distribution modelling of 180 of the 405 sources in the SIMBA source list. For each 
of the temperature, mass and luminosity, the range corresponding to a 1-sigma (68 per cent) probability of occurrence is presented, with 

the aiifl ,„,, r values roy)rosoi]tii],y; the lower aufl uy)por A'ahios of this laiigo. iosporti"\'oh'. 



Peak position Idcnt Fit Temperature Mass Luminosity 



RA 

(J2000) 




Dec 
{J2000) 




Source Name 

a, 6 


tracer 

c 


Type 
d 


T, 


K 


Tcoldmox 

K 




Mq 




Mmax 

M0 


Lq 




L0 


06 07 


46, 


,29 


-06 23 05.0 


G213.61-12.6 


mr 


SINGLE 


2, 


.8E+01 


3 


.lE+01 


9, 


,3E+02 


1, 


■4E+03 


7.6E+04 


3, 


■ 3E+05 


06 09 


06, 


,50 


+21 50 42.0 


G188.79+1.02 


r 


SED 


1, 


.2E+01 


3, 


.4E+01 


2, 


,0E+02 


1, 


■4E+03 


2.5E+02 


2, 


■ 9E+04 


09 03 


13, 


,50 


-48 


55 


30, 


,0 


G269. 45-1.47 


mr 


SED 


2, 


,7E+00 


3, 


.5E+01 


3, 


,lE+02 


7, 


■ 6E+03 


6.1E+02 


7, 


■ OE+04 


09 16 


41, 


,40 


-47 


55 


18, 


,0 


G270.25-I-0.84 


m 


SED 


2, 


.7E+00 


2, 


.9E+01 


1, 


.5E+02 


3, 


■8E+03 


8.9E+01 


1, 


■ lE+04 


10 48 


03, 


,98 


-58 


27 


04, 


,0 


G287.37-I-0.65 


m 


SED 


2, 


.7E+00 


3, 


.6E+01 


6, 


.6E+01 


2, 


■ 2E+03 


1.3E+02 


2, 


■OE+04 


10 57 


33, 


,00 


-62 


59 


10, 


,0 


G290.40-2.91 


m 


SED 


2, 


.7E+00 


2, 


.9E+01 


6, 


.6E+01 


1, 


■9E+03 


3.7E+01 


3, 


■9E+03 


11 12 


18, 


,06 


-58 


46 


19, 


,0 


G290.37-I-1.66 


m 


SED 


2, 


.7E+00 


3, 


,5E+01 


5, 


,0E+01 


1, 


■4E+03 


8.5E+00 


1, 


■ lE+04 


11 35 


31, 


,04 


-63 


14 


52, 


,3 


G294.52-1.6t 


m 


SED 


2, 


.7E+00 


3 


.2E+01 


1 


■ lE+Ol 


1, 


■ 5E+02 


l.lE+01 


1, 


■ lE+03 


12 11 


47, 


,64 


-61 


46 


29, 


,5 


G298.26+0.7 


m 


SED 


2, 


.7E+00 


3, 


,6E+01 


8, 


,7E+01 


2, 


■ 5E+03 


7.5E+01 


1, 


■ 7E+04 


12 35 


34, 


,95 


-63 


01 


35, 


,5 


G301. 14-0.2 


mr 


SED 


1 


.6E+01 


3, 


.9E+01 


1 


,lE+03 


3, 


■ lE+04 


6.0E+02 


2, 


■ 5E+05 


12 43 


32, 


,10 


-62 


55 


05, 


,8 


G302. 03-0.06 


mr 


SED 


2, 


.7E+00 


3, 


.6E+01 


2, 


,0E+02 


5, 


■ OE+03 


2.4E+02 


4, 


■ 5E+04 


13 10 


43, 


,25 


-62 


43 


04, 


,5 


G305.137-(-0.069t 


mm 


SINGLE 


3 


.4E+01 


3, 


.9E+01 


8, 


,7E+01 


1, 


■ 3E+02 


2.7E+04 


2, 


■ 6E+05 


13 16 


58, 


,38 


-62 


55 


25, 


,2 


G305.833-0.196t 


mm 


SINGLE 


2, 


.6E+01 


2, 


.8E+01 


7, 


.6E+01 


1, 


■OE+02 


4.2E+03 


1, 


■8E+05 


15 31 


44, 


,50 


-56 


30 


51, 


,0 


G323.74-0.30t 


m 


SED 


2, 


.7E+00 


3, 


.5E+01 


2, 


.7E+02 


6, 


■6E+03 


3.7E+02 


4, 


■8E+04 


16 11 


26, 


,90 


-51 


41 


57, 


,0 


G331.279-0.189 


m 


SED 


1, 


.7E+01 


3, 


.6E+01 


6 


.lE+02 


1, 


■3E+04 


8.2E+02 


1, 


■4E+05 


17 45 


54, 


,30 


-28 


44 


00, 


,0 


G0.204-t-0.05lt 


mm 


SINGLE" 


1, 


.3E+01 


3 


.5E+01 


4, 


.OE+02 


1, 


■9E+03 


5.2E+02 


2, 


■OE+05 


17 46 


04, 


,61 


-28 


24 


51, 


,0 


G0.49H-0.19t 


m 


SED" 


8, 


.lE+00 


3 


.2E+01 


5, 


.7E+01 


6, 


■ lE+02 


1.3E+01 


2, 


■ 7E+03 


17 46 


07, 


,09 


-28 


41 


28, 


,0 


G0.266-0.034t 


mm 


SINGLE" 


1, 


.2E+01 


3, 


.2E+01 


7, 


.lE+02 


3, 


■ 3E+03 


5.4E+02 


2, 


■ lE+05 


17 46 


07, 


,70 


-28 


45 


28, 


,0 


G0.21-0.00t 


mr 


SED" 


8, 


.4E+00 


3, 


.2E+01 


7 


.lE+02 


7, 


■ 6E+03 


9.8E+01 


3, 


■ 7E+04 


17 46 


08, 


,24 


-28 


25 


23, 


,0 


G0.497-l-0.170t 


mm 


SED" 


7, 


.OE+00 


2, 


.9E+01 


4, 


.3E+01 


5, 


■ 3E+02 


9.6E+00 


1, 


■ 8E+03 


17 46 


09, 


,52 


-28 


43 


36, 


,0 


G0.240-l-0.008t 


mm 


SINGLE" 


9, 


.4E+00 


2, 


.4E+01 


6, 


.6E+03 


3, 


■ lE+04 


1.8E+03 


2, 


■ 9E+05 


17 46 


10, 


,67 


-28 


23 


31, 


,0 


G0.527-|-0.18lt 


r 


SED" 


8, 


.lE+00 


3, 


.2E+01 


1, 


.lE+02 


1, 


■ 2E+03 


3.4E+01 


6, 


■6E+03 


17 46 


10, 


,74 


-28 


41 


36, 


,0 


G0.271-|-0.022t 


mm 


SINGLE" 


1, 


.OE+01 


3, 


.4E+01 


3, 


.lE+02 


1, 


■9E+03 


2.5E+02 


1, 


■8E+05 


17 46 


11, 


,35 


-28 


42 


40, 


,0 


G0.257-i-0.01lt 


mm 


SINGLE" 


9 


.8E+00 


2, 


.5E+01 


5, 


.7E+03 


2, 


■ 7E+04 


1.6E+03 


2, 


■9E+05 


17 46 


53, 


,41 


-28 


07 


27, 


,0 


G0.83-f0.18t 


m 


SED 


1, 


.3E+01 


4, 


.5E+01 


1, 


.OE+02 


9, 


■3E+02 


l.OE+02 


5, 


■3E+04 


17 47 


01, 


,19 


-28 


45 


36, 


,0 


G0.310-0.170t 


mm 


SINGLE" 


1, 


.6E+01 


6, 


.7E+01 


4, 


.3E+01 


2, 


■ 7E+02 


3.8E+02 


8, 


■ lE+05 


17 47 


09, 


,71 


-28 


46 


08, 


,0 


G0.32-0.20t 


mr 


SINGLE" 


9 


.8E+00 


2, 


.9E+01 


3 


.3E+03 


2, 


■OE+04 


1.4E+03 


4, 


■4E+05 


17 47 


20, 


,66 


-28 


46 


56, 


,0 


GO. 325-0.242+ 


mm 


SINGLE" 


1, 


.3E+01 


5, 


■ OE+Ol 


1 


.3E+02 


7, 


■ lE+02 


3.6E+02 


4, 


■4E+05 


17 48 


31, 


,59 


-28 


00 


30, 


,9 


Gl. 124-0.065 


mm 


SINGLE 


3, 


.4E+01 


5, 


.6E+01 


1, 


.7E+02 


4, 


■ OE+02 


1.2E+04 


7, 


■4E+05 


17 48 


34, 


,65 


-28 


00 


16, 


,0 


Gl. 134-0.073 


mm 


SINGLE" 


1, 


.2E+01 


5, 


.8E+01 


6, 


.6E+01 


5, 


■ 3E+02 


2.3E+02 


7, 


■ 2E+05 


17 48 


36, 


.41 


-28 


02 


31, 


,0 


Gl. 105-0.098 


mm 


SINGLE" 


8, 


.4E+00 


3, 


.OE+01 


9, 


.3E+02 


7, 


■6E+03 


6.3E+02 


4, 


■OE+05 


17 48 


42, 


,46 


-28 


01 


35, 


,0 


Gl.13-0.11 


r 


SED" 


3, 


.9E+00 


2, 


.OE+01 


2, 


.5E+03 


9, 


■3E+04 


2.5E+02 


2, 


■OE+05 


17 48 


49, 


,75 


-28 


01 


04, 


,0 


Gl.14-0.12 


m 


SINGLE" 


1, 


.5E+01 


6, 


.5E+01 


7, 


.6E+01 


5, 


■3E+02 


3.4E+02 


9, 


■3E+05 


17 50 


15, 


,11 


-27 


54 


23, 


,0 


GO.55-0.85 


mr 


SINGLE" 


9, 


.4E+00 


2, 


.9E+01 


6, 


.lE+02 


3, 


■8E+03 


2.5E+02 


1, 


■4E+05 


17 50 


18, 


,77 


-28 


53 


19, 


,0 


GO.549-0.868 


mm 


SINGLE" 


1, 


.4E+01 


6 


.2E+01 


7, 


.lE+00 


5, 


■OE+01 


4.9E+01 


2, 


■OE+05 


17 50 


25, 


,46 


-28 


50 


15, 


,0 


GO.627-0.848 


mm 


SINGLE" 


1, 


.5E+01 


7, 


.OE+01 


5, 


.3E+00 


3, 


■8E+01 


3.4E+01 


2, 


■OE+05 


17 50 


26, 


,07 


-28 


52 


31, 


,0 


GO.600-0.871 


mm 


SINGLE" 


1 


.OE+01 


4, 


.2E+01 


1, 


.9E+01 


1, 


■ 3E+02 


4.0E+01 


1, 


■4E+05 


17 50 


46, 


,50 


-26 


39 


44, 


,0 


G2.54+0.20t 


m 


SINGLE 


1 


.9E+01 


2 


.7E+01 


1, 


.7E+02 


4, 


■ OE+02 


5.6E+02 


1, 


■4E+05 


17 59 


02, 


,84 


-24 


20 


55, 


,0 


G5.48-0.24 


r 


SED" 


8 


.4E+00 


3, 


.4E+01 


1, 


.4E+03 


1, 


■ 5E+04 


5.3E+02 


1, 


■ 9E+05 


17 59 


07, 


,53 


-24 


19 


19, 


,0 


G5. 504-0.246 


mm 


SINGLE" 


9 


.8E+00 


3 


.2E+01 


8. 


.lE+02 


5, 


■ OE+03 


8.5E+02 


4, 


■ lE+05 


18 00 


30, 


,42 


-24 


03 


59, 


,0 


G5.89-0.39 


r 


SED" 


4, 


.lE+00 


1 


,9E+01 


7 


.lE+02 


2, 


■ 3E+04 


9.6E+01 


4, 


■ lE+04 


18 00 


40, 


,90 


-24 


04 


12, 


,0 


G5.90-0.42t 


m 


SED" 


9, 


.4E+00 


4, 


.OE+01 


3, 


.5E+02 


3, 


■8E+03 


1.3E+02 


6, 


■6E+04 


18 00 


49, 


,74 


-23 


20 


25, 


,0 


G6.53-0.10 


r 


SED 


1, 


.3E+01 


3, 


.5E+01 


3, 


.3E+03 


2, 


■ 7E+04 


5.6E+03 


6, 


■6E+05 


18 00 


54, 


,58 


-23 


16 


54, 


,0 


G6.60-0.08t 


m 


SED" 


4, 


.9E+00 


1, 


.2E+01 


1, 


.OE+00 


4, 


■OE+00 


4.9E-02 


4, 


■ 2E+01 


18 02 


49, 


,31 


-21 


48 


34, 


,0 


G8.111-I-0.257 


mm 


SINGLE" 


1, 


.OE+01 


6, 


.2E+01 


6, 


.lE+00 


7, 


■6E+01 


1.9E+01 


2, 


■ lE+05 


18 02 


52, 


,76 


-21 


47 


54, 


,0 


G8. 127+0.255 


mm 


SINGLE" 


5 


.9E+00 


3, 


.lE+01 


3, 


.8E+01 


6, 


■ lE+02 


2.2E+01 


1, 


■ 3E+05 


18 02 


56, 


,21 


-21 


47 


38, 


,0 


G8. 138+0.246 


mm 


SINGLE" 


6 


.8E+00 


3, 


.2E+01 


1, 


.OE+02 


1, 


■4E+03 


4.4E+01 


1, 


■ 3E+05 


18 03 


01, 


,95 


-21 


48 


02, 


,0 


G8. 13+0.22 


mr 


SED" 


6 


.5E+00 


3, 


■ OE+Ol 


5, 


.3E+02 


1, 


■ OE+04 


8.1E+01 


4, 


■ OE+04 


18 03 


26, 


,85 


-24 


22 


29, 


,0 


G5. 948-1. 125 


mm 


SINGLE" 


7 


.8E+00 


4, 


.7E+01 


3, 


.5E+00 


5, 


■ OE+01 


7.8E+00 


1, 


■ lE+05 


18 03 


29, 


,19 


-24 


21 


49, 


,0 


G5. 962-1. 128 


mm 


SINGLE" 


7 


.3E+00 


3, 


.2E+01 


7, 


■ lE+OO 


8, 


■ 7E+01 


9.3E+00 


8, 


■ 3E+04 


18 03 


33, 


,88 


-24 


21 


41, 


,0 


G5.975-1.146 


mm 


SINGLE" 


8, 


.7E+00 


5, 


.4E+01 


4, 


.6E+00 


5, 


■ 7E+01 


1.2E+01 


1, 


■ 5E+05 


18 03 


36, 


,80 


-24 


22 


08, 


,0 


G5.971-1.158 


mm 


SINGLE" 


1, 


.3E+01 


5, 


.4E+01 


1, 


.4E+01 


1, 


■OE+02 


7.0E+01 


2, 


■ lE+05 


18 03 


41, 


,49 


-24 


22 


37, 


,0 


G5.97-1.17 


r 


SED" 


3 


.2E+00 


1, 


.7E+01 


1, 


.3E+02 


5, 


■ 7E+03 


1.6E+01 


2, 


■OE+04 


18 05 


13, 


,33 


-18 


50 


30, 


,0 


G10.10+0.72t 


r 


SED 


5, 


.3E+00 


2, 


.OE+01 


1, 


.OE+00 


7, 


■ lE+00 


1.3E-01 


3, 


■3E+01 


18 06 


14, 


,80 


-20 


31 


29, 


,0 


G9. 63+0. 19 


mr 


SED" 


7 


.8E+00 


3, 


.lE+01 


2, 


.7E+02 


2, 


■ 8E+03 


5.2E+01 


1, 


■ 3E+04 


18 06 


18, 


,91 


-21 


37 


21, 


,0 


G8. 68-0.36 


mr 


SINGLE" 


8 


.4E+00 


2, 


.5E+01 


3, 


.3E+03 


2, 


■ OE+04 


7.3E+02 


2, 


■ 2E+05 


18 06 


23, 


,49 


-21 


36 


57, 


,0 


G8. 686-0. 366 


m 


SINGLE" 


1 


.lE+01 


3, 


.lE+01 


7, 


.lE+02 


3, 


■ 8E+03 


4.7E+02 


2, 


■ OE+05 


18 07 


50, 


,36 


-20 


18 


51, 


,0 


G9. 99-0.03 


m 


SINGLE" 


8 


.7E+00 


2, 


.7E+01 


3, 


.lE+02 


1, 


■ 9E+03 


9.8E+01 


9, 


■ 7E+04 


18 07 


53, 


,21 


-20 


18 


19, 


,0 


GlO.001-0.033 


r 


SED" 


7 


.OE+00 


3, 


.4E+01 


5, 


■ OE+Ol 


8, 


■ lE+02 


9.9E+00 


2, 


■ 5E+03 


18 08 


38, 


,47 


-19 


51 


48, 


,0 


G10.47+0.02t 


mr 


SED" 


7, 


.3E+00 


2, 


.6E+01 


7, 


.6E+03 


9, 


■3E+04 


7.1E+02 


1, 


■ 5E+05 


18 08 


45, 


,47 


-19 


54 


30, 


,0 


G10.44^0.0lt 


m 


SINGLE" 


1, 


.3E+01 


4, 


.2E+01 


3 


.5E+02 


1, 


■9E+03 


6.4E+02 


4, 


■3E+05 
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Table 2 - continued 



Peak position Ident Fit Temperature Mass Luminosity 



RA 

(J2000) 




Dec 
(J2000) 




Source Name 
o,6 


tracer 

c 


Type 

d 


T, 


Zoldmin. 

K 


Tcoldmax 

K 




Mq 




M© 




Emin 

L© 




^max 

L© 


18 08 


45. 


,85 


-20 


05 


42. 


,0 


GlO.287-0.110' 


mm 


SINGLE 


7. 


.8E-I-00 


2. 


.6E-I-01 


8. 


,7E-|-01 


6 


.lE+02 


6. 


,1E+01 


9 


.4E+04 


18 08 


49, 


.25 


-20 


05 


58, 


.0 


GIO. 284-0. 126 


m 


SED™ 


7, 


.8E-I-00 


3 


.4E-I-01 


5. 


.7E-I-0I 


7, 


.IE+02 


2. 


OE+OI 


7 


.lE+03 


18 08 


52, 


.66 


-20 


05 


58, 


.0 


GIO. 288-0. 127' 


mm 


biNGLii," 


1 


.3E-I-01 


4, 


.8E-I-01 


1. 


.6E+01 


1 


.OE+02 


6. 


.2E+01 


1 


.4E+05 


18 08 


56, 


.07 


-20 


05 


50, 


.0 


GIO. 29-0. 14 


mr 




7, 


.3E-I-00 


2, 


.9E-I-01 


2. 


.OE+02 


2 


.5E+03 


5, 


.2E+01 


1 


.5E+04 


18 09 


00, 


.04 


-20 


03 


34, 


.0 


CjiU.o4o-U.i4zi 


m 


biJNGLii/ 


1 


.3E-I-01 


4, 


.2E-I-01 


4. 


.3E+01 


2 


.OE+02 


8, 


.2E+01 


1 


. lE+05 


18 09 


03, 


.49 


-20 


02 


54, 


.0 


GIO. 359-0. 149' 


mm 


blNGLlj 


1 


. lE-fOl 


4, 


.7E-I-01 


2. 


.2E+01 


1 


.7E+02 


3, 


.7E+01 


1 


.2E+05 


18 09 


21, 


.03 


-20 


19 


25, 


.0 


GIO. 15-0. 34 


r 


SINGLE 


3, 


.8E+00 


1, 


.9E+01 


9. 


.3E+02 


2 


.3E+04 


1, 


.8E+02 


4 


.lE+05 


18 10 


15, 


.59 


-19 


54 


45, 


.0 


G10.63-0.33B 


mm 


SINGLE'* 


1 


.lE-l-01 


5, 


.2E-I-01 


2. 


.OE-l-02 


1, 


.9E+03 


5, 


.5E+02 


1 


.lE+06 


18 10 


18, 


.42 


-19 


54 


29, 


.0 


GIO. 62-0.33 


m 


SED™ 


7, 


.OE-t-00 


3, 


.OE-hOl 


8. 


.lE-l-02 


1, 


.lE+04 


2, 


.OE+02 


4, 


.6E+04 


18 10 


19, 


.00 


-20 


45 


25, 


.0 


G9.88-0.75 


r 


SED™ 


7, 


.8E-I-00 


2, 


.6E-I-01 


8. 


. lE-l-02 


7, 


.6E+03 


8, 


.9E+01 


1 


.6E+04 


18 10 


23, 


.56 


-20 


43 


09, 


.0 


G9. 924-0. 749 


mm 


SED™ 


4, 


.9E-I-00 


2, 


.5E-I-01 


5. 


.7E-I-01 


1 


.lE+03 


1, 


.6E+01 


7 


.3E+03 


18 10 


28, 


.77 


-19 


55 


48, 


.0 


GIO. 62-0. 38 


mr 


biNGLE*^ 


8, 


.7E-I-00 


2, 


.8E+01 


8. 


.7E+03 


6 


.lE+04 


2, 


.3E+03 


8 


.OE+05 


18 11 


23, 


.87 


-19 


32 


20, 


.0 


Gil. 075-0. 384 


mm 


SED" 


4, 


.9E-I-00 


2, 


.5E-f01 


5. 


.7E+01 


1 


.lE+03 


1, 


.6E+01 


7. 


.3E+03 


18 11 


31, 


.80 


-19 


30 


44, 


.0 


/~1 -1-1 -1-1 r\ A 

Gil. 11-0.34 


r 


SED" 


8, 


.lE-l-00 


2, 


.8E-I-01 


8, 


.lE+02 


7, 


.6E+03 


1, 


.2E+02 


2 


.8E+04 


18 11 


35, 


.76 


-19 


30 


44, 


.0 


Cjil.il7-U.41o 


mm 


CTMi^T TT'CK 

biJNGLlif 


1 


.2E-t-01 


3, 


.6E-I-01 


1, 


.3E+02 


6 


. lE+02 


1, 


.3E+02 


1 


.2E+05 


18 11 


51, 


.40 


-17 


31 


30, 


.0 


G12.88+0.48' 


m 


SIN GLE 


8, 


.lE-l-00 


2, 


.6E+01 


1, 


.lE+03 


7, 


.6E+03 


2, 


.6E+02 


1 


.4E+05 


18 11 


52, 


.90 


-18 


36 


03, 


.0 


Gil. 948-0. 003' 


mm 


SED 


2, 


.7E-I-00 


4, 


.3E-I-01 


1, 


.5E-I-02 


3, 


.8E+03 


6, 


.9E+02 


1 


.2E+05 


18 11 


53, 


.64 


-17 


30 


02, 


.0 


G12. 914-1-0. 493 ' 


mm 


SIN GLE" 


6 


C5TT^ 1 t\f\ 

.8EH-00 


2, 


.9E-I-01 


6, 


.6E-f01 


7, 


.lE+02 


3, 


.5E+01 


9 


.OE+04 


18 12 


01, 


.34 


-18 


31 


55, 


.0 


Gi2.O2-O.O3t 


m 


SED 


1, 


.4E-I-01 


2, 


.7E-I-01 


2, 


.3E-I-02 


9, 


.3E+02 


2, 


.2E+02 


1 


. lE+04 


18 12 


11, 


.11 


-18 


41 


30, 


.0 


Gil. 903-0. 140' 


mr 


SIN GLE" 


4, 


.7E-I-00 


1, 


.7E-I-01 


3, 


.5E-I-02 


5 


.7E+03 


3, 


.3E+01 


6 


.5E+04 


18 12 


15, 


.61 


-18 


44 


58, 


.0 


Gil. 861-0. 183^ 


mm 


biNGLE" 


1 


.4E+01 


5 


.6E+01 


1, 


.2E+01 


7, 


.6E+01 


5, 


.3E+01 


1 


.5E+05 


18 12 


17, 


.30 


-18 


40 


02, 


.0 


Cjil.y0-U.i4i 


m 


biiNGLE 


I 


.OE+01 


2, 


.6E+01 


1, 


.5E+02 


7, 


. lE+02 


8, 


.5E+01 


8 


.5E+04 


18 12 


19, 


.55 


-18 


39 


54, 


.0 


Gil. 942-0. 157' 


mm 


SINGLE 


8 


.7E-I-00 


3, 


AE+01 


7, 


.6E+01 


7, 


. lE+02 


5, 


.OE+01 


9 


.6E+04 


18 12 


33, 


.13 


-18 


30 


05, 


.0 


10 110 r\ lor^ 

G12. 112-0. 125 


mm 


biNGLE 


2 


.5E+01 


2, 


.7E4-01 


1, 


.2E+03 


1 


.6E+03 


4, 


.lE+04 


2 


.8E+05 


18 12 


39, 


.31 


-18 


24 


13, 


.0 


Gi2.20-0.09 


mr 


biJNGLE 


1 


.OE-f 01 


3, 


.0E-(-01 


8, 


.7E-f 03 


4 


.6E+04 


4, 


.2E+03 


9 


.2E+05 


18 12 


43, 


.25 


-18 


25 


09, 


.0 


10 10 r\ 10A 

G12.18-0.12A 


m 


SIN GLE 


1 


.3E-I-01 


4, 


.3E-I-01 


8, 


.lE-l-02 


5 


.OE+03 


1, 


.3E+03 


8, 


.9E+05 


18 12 


44, 


.37 


-18 


24 


21, 


.0 


G12. 216-0. 119 


mm 


<!^T1VT/^T TTl CV 

SINGLE" 


1 


.2E-I-01 


3, 


.6E-I-01 


1, 


.9E-I-03 


1, 


.OE+04 


2, 


.6E+03 


8, 


.4E+05 


18 12 


50, 


.64 


-18 


40 


31, 


.0 


Gil. 99-0.27' 


m 


SED™ 


8, 


.7E-I-00 


3, 


.2E-I-01 


7, 


.6E-I-01 


7, 


.lE+02 


1, 


.5E+01 


3 


.6E+03 


18 12 


54, 


.72 


-18 


11 


04, 


.0 


G12.43-0.05 


r 


SED 


1 


.lE-l-01 


3 


.OE-l-01 


1, 


.9E-I-03 


1 


.3E+04 


1, 


.lE+03 


1 


.6E+05 


18 13 


54, 


.14 


-18 


01 


41, 


.0 


G12.68-0.18t 


m 


SED" 


7. 


.3E+00 


2 


.8E+01 


9, 


.3E+02 


8 


.7E+03 


2, 


.5E+02 


2 


.8E+04 


18 13 


58, 


.08 


-18 


54 


14, 


.0 


Gil. 94-0. 62B 


mm 


SINGLE 


6 


.8E+00 


1, 


.9E+01 


7, 


.lE+02 


5 


.OE+03 


1, 


.6E+02 


1 


.OE+05 


18 14 


00, 


.90 


-18 


53 


18, 


.0 


Gil. 93-0. 61 


mr 


SED™ 


7 


.8E-f00 


3, 


.OE-l-01 


6, 


.lE+02 


7, 


.6E+03 


7, 


.4E+01 


2 


.8E+04 


18 14 


07, 


.04 


-18 


00 


37, 


.0 


in '~7 n 0101 

G12. 722-0. 218 ' 


mm 


SED™ 


5 


.9E+00 


2, 


.lE-l-01 


3, 


.lE+02 


4, 


.3E+03 


9, 


.7E+01 


1 


.6E+04 


18 14 


33, 


.90 


-17 


51 


44, 


.0 


10 nn r\ or* 'P* T 

G12.90-0.25B ' 


mm 


SINGLE 


8 


.4E+00 


3, 


.6E+01 


1, 


.3E+02 


1 


.2E+03 


6, 


.lE+01 


1 


.3E+05 


18 14 


35, 


.54 


-16 


45 


36, 


.0 


10 0^1 n 

G13. 87-1-0. 28 


m 


SED 


9, 


.lE-l-00 


3, 


.7E-I-01 


5, 


.3E-I-02 


7, 


.6E+03 


3, 


.OE+02 


1 


.2E+05 


18 14 


36, 


.13 


-17 


54 


56, 


.0 


G12. 859-0. 272t 


mm 


SED™ 


6 


.5E-I-00 


3, 


.2E-I-01 


1, 


.7E-I-02 


3, 


.3E+03 


3, 


.4E+01 


1 


.OE+04 


18 14 


38, 


.94 


-17 


51 


52, 


.0 


G12.90-0.26t 


m 


SED™ 


7, 


.3E-t-00 


2, 


.9E-I-01 


9, 


.3E-I-02 


1, 


.3E+04 


1, 


.3E+02 


3 


.3E+04 


18 16 


22, 


.10 


-19 


41 


19, 


.0 


Gil. 49-1. 48' 


m 


SED 


1 


.lE-f-Ol 


2, 


.8E-I-01 


4, 


.3E-f01 


3, 


.lE+02 


1, 


.3E+01 


2 


.2E+03 


18 17 


02, 


.17 


-16 


14 


28, 


.0 


G14.dO-I-0.01 ' 


mr 


SED™ 


8, 


.lE-l-00 


2, 


.9E-I-01 


1, 


.5E-I-02 


1, 


.6E+03 


2, 


lE+Ol 


4, 


.3E+03 


18 19 


12, 


.03 


-20 


47 


23, 


.0 


GIO. 84-2. 59 


r 


OTAT/~1T TTl rv 

SIN GLE" 


7, 


.8E-I-00 


2, 


.5E-I-01 


1, 


.5E-I-02 


1 


.IE+03 


6, 


.6E+0I 


8 


.3E+04 


18 20 


23, 


.10 


-16 


11 


31, 


.0 


G15.03-0.67' 


mr 


SIN GLE 


8 


.4E+00 


3 


.4E+01 


1, 


.lE+03 


1 


.OE+04 


4, 


.7E+02 


4 


.OE+05 


18 21 


09, 


.10 


-14 


31 


40, 


.0 


G16.58-0.05' 


m 


SED™ 


7 


.6E+00 


2, 


.6E+01 


6, 


.lE+02 


6 


.6E+03 


3, 


.4E+01 


9 


.4E+03 


18 21 


14, 


.61 


-14 


32 


52, 


.0 


G16. 580-0. 079' 


mm 


SIN GLE 


5 


.lEH-00 


1, 


.OE-l-01 


4, 


.OE+02 


2 


.2E+03 


1, 


.2E+01 


3 


.2E+04 


18 25 


01, 


.30 


-13 


15 


27, 


.0 


G18. 15-0.28 


r 


SIN GLE 


3 


.2E-I-00 


1, 


.7E-I-01 


5, 


.7E+01 


2 


.2E+03 


1, 


.lE+01 


8, 


.7E+04 


18 25 


07, 


.33 


-13 


14 


23, 


.0 


G18. 177-0.296 


mm 


OTTVT/~1T TTl rv 

SIN GLE" 


1 


.OE-l-01 


4, 


.8E-I-01 


2, 


.8E+01 


2, 


.7E+02 


4, 


.7E+01 


1 


.6E+05 


18 25 


41, 


.65 


-13 


10 


16, 


.0 


G18. 30-0.39 


r 


SED 


9, 


.8E-t-00 


3, 


.5E-I-01 


2, 


.OE+02 


3 


.3E+03 


6, 


.5E+01 


3 


.3E+04 


18 27 


16, 


.34 


-11 


53 


51, 


.0 


G19.61-0.lt 


m 


SED™ 


8, 


.4E-t-00 


3, 


.4E-I-01 


1, 


.5E+02 


1, 


.9E+03 


1, 


.4E+01 


7, 


.8E+03 


18 27 


37, 


.86 


-11 


56 


40, 


.0 


G19. 607-0. 2341 


mr 


SED 


7, 


.8E-I-00 


3, 


.5E-I-01 


5, 


.3E+02 


1, 


.3E+04 


1, 


.7E+02 


7, 


.9E+04 


18 27 


55, 


.30 


-11 


52 


48, 


.0 


G19.70-0.27Ai 


m 


SED 


1 


.7E-I-01 


3 


.9E-I-01 


5, 


.3E+02 


1 


.3E+04 


1, 


.6E+03 


2 


.3E+05 


18 29 


24, 


.20 


-15 


16 


06, 


.0 


G16.86-2.15' 


m 


SED™ 


6 


.8E-I-00 


1, 


.8E-I-01 


7, 


.lE+02 


5 


.OE+03 


5, 


.6E+01 


2 


.6E+03 


18 29 


33, 


.60 


-15 


15 


50, 


.0 


G16.883-2.188t 


mm 


SINGLE™ 


1 


.OE-l-01 


3, 


.2E-f01 


1, 


.lE+01 


6, 


.6E+01 


9, 


.4E+00 


4 


.OE+04 


18 31 


02, 


.64 


-09 


49 


38, 


.0 


G21.87-l-0.0lt 


mr 


SED™ 


8 


.4E-I-00 


3, 


.2E+01 


2, 


.8E+01 


2, 


.7E+02 


8, 


.8E+00 


2 


.lE+03 


18 31 


43, 


.02 


-09 


22 


28, 


.0 


G22.36-I-0.07 


m 


SED™ 


2, 


.7E-I-00 


1, 


.3E+01 


1, 


.7E+02 


1, 


.3E+04 


1, 


.4E+01 


3 


.3E+03 


18 33 


53, 


.06 


-08 


07 


23, 


.0 


G23.71-I-0.17 


r 


SED™ 


7, 


.3E-I-00 


2, 


.8E-I-01 


1, 


.lE+03 


1, 


.lE+04 


1, 


.7E+02 


5 


.9E+04 


18 33 


53, 


.60 


-08 


08 


51, 


.0 


G23.689-I-0.159 


mm 


SINGLE™ 


1, 


.2E-I-01 


4, 


.7E-I-01 


8, 


.7E+0I 


4, 


.6E+02 


2, 


.lE+02 


2 


.9E+05 


18 34 


09, 


.23 


-07 


17 


45, 


.0 


G24.47-l-O.49t 


r 


SED™ 


9 


.lEH-00 


3, 


.7E-I-01 


8, 


.lE+02 


1, 


.OE+04 


1, 


.9E+02 


1 


.9E+05 


18 34 


20, 


.90 


-05 


59 


48, 


.0 


G25.65-I-1.04 


mr 


SED™ 


7, 


.8E-t-00 


2, 


.6E-I-01 


6, 


.lE+02 


5, 


.7E+03 


8, 


.lE+01 


1 


.3E+04 


18 34 


31, 


.30 


-08 


42 


47, 


.0 


G23.25-0.24t 


m 


SINGLE" 


1, 


.OE-l-01 


2, 


.8E-I-01 


7, 


.6E+01 


3 


.5E+02 


3, 


.5E+01 


6 


.6E+04 


18 34 


36, 


.16 


-08 


42 


39, 


.0 


G23.268-0.257t 


mm 


SINGLE™ 


5 


.4E+00 


1, 


.7E+01 


8, 


.lE+02 


8, 


.7E+03 


5, 


.4E+01 


7. 


.2E+04 


18 34 


39, 


.20 


-08 


31 


41, 


.0 


G23.43-0.18t 


m 


SINGLE™ 


1 


.OE+01 


2, 


.5E+01 


1, 


.9E+03 


8, 


.7E+03 


7, 


.6E+02 


1 


.8E+05 


18 36 


06, 


.69 


-07 


13 


47, 


.0 


G23.754+0.095t 


mm 


SINGLE™ 


9 


.lE+00 


3, 


.2E+01 


3, 


.lE+02 


2 


.5E+03 


1, 


.6E+02 


1 


.5E+05 


18 36 


12, 


.60 


-07 


12 


11, 


.0 


G24.78-l-0.08t 


m 


SINGLE™ 


7, 


.8E-I-00 


2, 


.3E-I-01 


5, 


.7E+03 


4, 


.OE+04 


9, 


.5E+02 


2 


.4E+05 


18 36 


17, 


.86 


-07 


08 


52, 


.0 


G24.84-f0.08t 


m 


SINGLE" 


1, 


.OE-fOl 


3 


.4EH-01 


5, 


.3E+02 


3 


.3E+03 


3, 


.6E+02 


2 


.OE+05 



8 T. 

Table 2 



Hill et al. 

continued 



Peak position Ident Fit Temperature Mass Luminosity 



RA 

(J2000) 




Dec 
(J2000) 


Source Name 

a.h 


tracer 

c 


Type 

d 


T< 


K 


Tcoldn 


" K 




Mq 




Mq 




L© 




L© 


18 


36 


25, 


,92 


-07 05 16.9 


G24, 


,919+0.088t 


mm 


SED 


9, 


,8E+00 


3, 


,2EH 


-01 


8 


,lE+02 


8 


,7E+03 


2, 


,9E+02 


8 


,5E+04 


18 


38 


03, 


,00 


-06 24 09.0 


G25, 


,70-f0.04 


mr 


SED" 


7, 


,0E+00 


3, 


,OEH 


-01 


1 


,9E+03 


2 


,7E+04 


2, 


,3E+02 


7, 


,9E+04 


18 


39 


03, 


,94 


-06 24 13.0 


G25, 


,82-0.17t 


m 


blNGLhj" 


8, 


,7E+00 


2, 


,3EH 


-01 


1 


,9E+03 


1 


,0E+04 


4, 


,4E+02 


1 


,5E+05 


18 


42 


42, 


,60 


-04 15 39.0 


G28, 


,14-O.OOt 


m 


SED 


1, 


,3E+01 


3, 


,7EH 


-01 


2 


,0E+02 


1 


,4E+03 


1, 


,0E+02 


3 


,9E+04 


18 


42 


54, 


,89 


-04 07 40.0 


G28, 


,287-fO.OlOt 


mm 




1, 


,0E+01 


3, 


,9EH 


-01 


1 


,0E+02 


8 


, lE+02 


1, 


,4E+02 


2 


,lE+05 


18 


42 


58, 


,10 


-04 13 56.0 


G28, 


,20-0.04t 


mr 


SED" 


7, 


,8E+00 


2, 


,8EH 


-01 


2 


,2E+03 


2 


,7E+04 


2, 


,5E+02 


6 


,3E+04 


18 


43 


02, 


,91 


-04 14 52.0 


G29, 


,193-0.073t 


mm 


SINGLE 


1, 


,0E+01 


4, 


,3EH 


-01 


6 


,6E+01 


5 


,3E+02 


1, 


,3E+02 


2 


,4E+05 


18 


44 


15, 


,17 


-04 01 56.0 


G28, 


,28-0.35 


mr 


SED" 


5, 


,7E+00 


2, 


,9EH 


-01 


1 


,7E+02 


4, 


,3E+03 


3, 


,2E+01 


2 


,2E+04 


18 


44 


21, 


,57 


-04 17 35.0 


G28, 


,31-0.38t 


m 


SED" 


9, 


,4E+00 


4, 


,3EH 


-01 


1 


,3E+02 


1 


,6E+03 


8, 


,2E+01 


7, 


,lE+04 


18 


45 


52, 


,76 


-02 42 29.0 


G29, 


,888-fO.OOlt 


mm 


SINGLE" 


5, 


,4E+00 


2, 


,3EH 


-01 


2 


,7E+02 


4, 


,3E+03 


1, 


,lE+02 


3 


,3E+05 


18 


45 


54, 


,36 


-02 42 45.0 


G29, 


,889-0.006t 


mm 


SINGLE 


5, 


,lE+00 


1, 


,6EH 


-01 


2 


,3E+02 


2 


,8E+03 


3, 


,0E+01 


1 


,2E+05 


18 


45 


59, 


,70 


-02 41 17.0 


G29, 


,918-0.014t 


mm 


SINGLE" 


6, 


,8E+00 


3, 


,0E+01 


6 


,6E+01 


8 


,lE+02 


7, 


,7E+01 


3 


,8E+05 


18 


46 


00, 


,23 


r\i^ A ^ r\r\ r\ 

-02 45 09.0 


G29, 


,86-0.04t 


m 


SED" 


1 


,0E+01 


4, 


,2EH 


-01 


2 


,7E+02 


2 


,5E+03 


1, 


,lE+02 


4, 


,3E+04 


18 


46 


01, 


,30 


-02 45 25.0 


G29, 


,861-0.053t 


mm 


SINGLE 


1 


,0E+01 


5, 


,6EH 


-01 


7, 


,6E+01 


8 


, lE+02 


2, 


,8E+02 


1 


,0E+06 


18 


46 


02, 


,37 


-02 45 57.0 


G29, 


,853-0.062t 


mm 


SINGLE" 


1 


,2E+01 


5, 


,4EH 


-01 


1 


,3E+02 


1 


,lE+03 


5, 


,0E+02 


1 


,0E+06 


18 


46 


03, 


,97 


nr\ or" n 

-02 39 25.0 


G29, 


,96-0.02Bt 


mr 


SED" 


7, 


,3E+00 


3, 


,4EH 


-01 


1 


,9E+03 


3 


, lE+04 
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" Distance values arc reported in |Hill eiTaT} |2005| )- A in thits column indicates those sources which have a distance ambiguity. 

^ Source names given to two (or less) decimal places are consistent with those reported by lWalsh et al.| U^^d l ; fCTinier et all jj^OOU ) ; fThompson eT'm (^HH^) which were 
targeted in the SIMBA survey (Paper I). Note that in a few instances, there may not be a direct translation between their Galactic name and their equatorial 
coordinates. Source names given to three decimal places, denote those sources identified by SIMBA, with the extended Galactic names intended to distinguish closely 
associated sources. The source names are consistent with Table 5, Paper I. 

^ Denotes (any) association with methanol maser sites (m) and/or UC Hll regions (r). Those sources associated with both a methanol maser and a radio continuum 
source are depicted by (mr). The MM-only sources which are devoid of both of these sources are denoted (mm). 

This column depicts the type of SED applied to the source. 'SED' indicates that a two-component spectral energy distribution has been applied, whilst 'SINGLE' 
indicates that a single cold-component fit has been applied, to the source. An in this column denotes those sources that were fit with IRAS upper limits. 
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Figure 1. SED fits and probability distributions of the temperature and the mass for three sources, from top to bottom: G 12.02 - 0.03, 
GO. 21-0. 00 and G 10.288 - 0.127. Left: example spectral energy distributions. The blue and red lines represent the hot and cold compo- 
nents, respectively. The black line represents the total SED fit. Crosses are observational data points whilst triangles are upper limits. 
Centre and right: histograms of marginal probability distributions for the temperature and mass of the cold component, respectively. 
The shaded component represents the range of validity of the parameters presented in Table [2] This is the range of values that encloses 
the 68 per cent probability (i.e. the shaded area represents 0.68 of the area under the curve) as defined in Equation [2] 



ambiguity. The results were consistent with each other and 
indicated the small influence of assuming the near distance 
value for those with an ambiguity. We thus refrain from pre- 
senting the far distance determinations of each of the pa- 
rameters in this work and refer the reader to Paper I for the 
far distances, and Equation[l]for the appropriate scaling fac- 
tor where the far distance value is more applicable. Sources 
with a distance ambiguity are denoted by a Mn Table (2] 

Table [5] presents the temperature, mass and luminosity 
of each of the 180 sources for which SED models were drawn. 
Figure [1] provides an example SED and probability curves 
for both the temperature and mass of three sources. The top 
and middle panels of this figure display the respective plots 
for a two component SED, whilst the bottom panel illus- 
trates a single component fit. The top panel presents a well 
constrained source with sharply peaked probability curves 
for both the temperature and mass. The middle and bottom 
panels on the other hand present sources with less well con- 
strained data, i.e. the SED has had the IRAS upper limits 
incorporated. Figure [1] illustrates the need for observations 



spanning a broad range of the wavelength parameter space of 
the SED, as well as ample and tightly constrained data for 
obtaining sharp probability curves, i.e. strong constraints 
on the parameters (for instance G 12.02 -0.03, top panel). 
When less data are available (for instance, G 0.21 -0.00 and 
G 10.288-0.127, middle and bottom panels, respectively) for 
SED analysis, the resultant probability curves are wider and 
uncertainties on the derived parameters are larger. In this 
instance, estimating the model parameters from the best fit 
values is simply not accurate enough. 

Class-comparative histogram and cumulative distribu- 
tion plots for the temperature, mass and luminosity of the 
sample are usually built by sorting the individual sources 
into bins according to the values for parameters obtained for 
the best SED models. However, we have shown that these 
"best values" can be attached to significant uncertainties 
that should be taken into account when doing statistical 
studies of the complete sample. 

For a given source, instead of incrementing the parame- 
ter bins corresponding to the "best value" by 1, we increment 
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all bins by the Bayesian probability that the parameter takes 
this value, with the sum of all probabilities being 1. For well 
constrained sources, the probability is sharp and peaks at 
the "best value" (see. Fig.[l|. In this instance our procedure 
is almost equivalent to the classical one. For less well con- 
strained sources however, a larger number of bins are incre- 
mented by a small amount. With this method, uncertainties 
on the derived parameters of individual sources are auto- 
matically taken into account in the cumulative distribution 
plots and more weight is given to well constrained sources. 
That is, less well constrained sources essentially have fiat 
distributions and they do not contribute to the shape of the 
histogram or cumulative plots. 

Cumulative distribution plots for each parameter are 
presented in Fig. (5] For comparative purposes, these distri- 
butions depict each of the four classes of source (see sec- 
tion [27l]) on a single plot for each parameter. In an attempt 
to ascertain whether the MM-only cores have characteris- 
tics similar to sources with star formation activity, we also 
drew cumulative distributions comparing the MM-only sam- 
ple with the star formation activity sample (the combination 
of classes M, MR and R sources) for each parameter. 

As we are primarily interested in the dust properties, 
temperature comparisons amongst the sources are made us- 
ing the Tcoid component of the sources, rather than the T^ot 
value which has only been determined for those sources with 
a full two-component SED. The mean and median values of 
the temperature, mass and luminosity for each of the classes 
of source, as well as the combined maser and radio contin- 
uum sources (class M, MR and R) and the entire sample are 
presented in Table [l] 

In order to ascertain the influence of the different types 
of fit - single or two-component - on the cumulative distribu- 
tions, we separated these two populations and drew cumu- 
lative distributions for each parameter for each fit type. The 
cumulative plots of each component (i.e. fit type) displayed 
little difference from the composite data presented in Fig. [2] 
for each of the temperature, mass and luminosity, and are 
thus not presented here. 

Kolmogorov-Smirnov (KS) tests were also performed 
from the cumulative distributions, in order to test the hy- 
pothesis that the MM-only sources are drawn from the same 
parent distribution as the other three classes of source (M, 
MR and R) . KS tests were applied to the temperature, mass 
and luminosity. The results from the KS tests can be found 
in Table[5] with discussion of each parameter in the following 
subsections. 

Generally, in order to conclude that two distributions 
are not drawn from the same sample, the KS probabilities 
must be small, 0.01. As a calibration measure, the KS 
tests were applied to class M and R sources for each of the 
parameters mentioned above. The results indicate that the 
likelihood of these two classes of source being from different 
distributions, for all three parameters tested, is small (see 
Table E}. 

4.2 Temperature 

The cumulative distribution plot of the temperature (Fig. [S] 
top, left) indicates that there is little distinction between the 
four classes of source in the sample. The sources with radio 
continuum associations (classes MR and R) appear to be the 



Table 5. Results from the KS test of temperature, mass and lu- 
minosity for the four classes of source as well as the star formation 
activity sample. Column 1 indicates the parameter being tested. 
Columns 2 and 3 list the classes of source being tested. Column 
4 gives the resultant KS probability that the objects in columns 
2 and 3 are from the same parent distribution. If this probability 
is <0.01 it is generally concluded that the samples are not drawn 
from the same population. 
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coolest of the sample, whilst conversely the MM-only cores 
are the warmest of the sample. We caution however that this 
distinction is very slight, as indicated by the KS tests which 
confirm that these data do not allow us to strongly discrim- 
inate between the different classes of source with respect to 
their temperature. 

The median temperature of the sample (i.e. the value 
corresponding to the 0.5 fraction value on Fig. [2]) of the MM- 
only sample is 17 K, the maser sample (class M) is 15 K, the 
maser+radio sample (class MR) is 13 K, and the radio sam- 
ple (class R) is 13 K. These median values are not incon- 
sistent with the 20 K temperature assumed in Paper I for 
purposes of mass derivation and they emphasize the small 
difference between the samples in terms of temperature. The 
cumulative distribution of the sources with known star for- 
mation activity (classes M, MR and R) indicates that collec- 
tively these sources are marginally cooler on average than 
class MM-only sources i.e. the sources apparently without 
star formation activity. 

The shape of each of the individual class distributions in 
a cumulative plot is also a useful diagnostic tool. The linear 
shape of the temperature profile of classes MR and R sources 
indicates that there is little constraint on their temperature 
as determined from our SED analysis. Only an upper limit of 
40 - 50 K can be firmly established. We also draw attention to 
the fact that the cumulative distribution of the temperature 
shows values <5K. As the cumulative distribution reports 
the sum of all probabilities over all temperatures, this is 
simply the tail end of the temperature distribution and we 
do not attribute any significance to these low temperatures. 

Interestingly, the cumulative distribution of the maser 
sample (class M) directly traces the temperature profile of 
the entire sample (not shown). This perhaps indicates that 
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Figure 2. Cumulative distributions of the SIMBA sources. The distributions of the individual classes are displayed on each of the plots, 
with the key on the top-left plot indicating which class of source is represented by which distribution, top left : distribution of temperature. 
top right: distribution of mass, bottom left: 2-dimensional probability distribution of the temperature and the mass parameters, bottom 
right: distribution of luminosity. 



the maser sample displays a global temperature profile typi- 
cal of a large cross section of massive star formation sources, 
and more specifically our entire sample. That is, the maser 
population may be considered the 'standard' massive star 
formation population (this is explored further in section [F!3)) . 



4.3 Mass 

The cumulative distribution plot of the mass (Fig. [2j top, 
right) indicates that MM-only sources are the least mas- 
sive of the sample, whilst the sources with both a methanol 
maser and radio continuum source (class MR) are the most 
massive of the sample. As per the temperature distribution. 
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Table 4. Mean and median as derived from the SED analysis of the four classes of source. For a breakdown of the number of sources 
within each class for which different SEDs were applied, refer to Table |3]. 











Class 




whole 


all except MM-only 


Parameter 




MM-only(MM) 


maser(M) 


maser+radio(MR) 


radio(R) 


sample 


(M+MR+R) 


temperature (K) 


mean 


16.0 


13.8 


12.5 


12.4 


14.2 


13.1 




median 


16.6 


14.9 


13.2 


13.0 


15.0 


13.9 


mass (M0) 


mean 


4.6E+02 


9.6E+02 


3.3E+03 


2.0E+03 


9.5E+02 


1.6E+03 




median 


4.5E+02 


8.9E+02 


3.0E+03 


2.2E+03 


9.4E+02 


1.6E+03 


luminosity (Lq) 


mean 


5.6E+03 


2.9E+03 


4.8E+03 


3.3E+03 


4.2E+03 


3.4E+03 




median 


5.2E+03 


3.3E+03 


5.7E+03 


5.4E+03 


4.7E+03 


4.4E+03 



the maser sample again traces the mass distribution of the 
whole sample (not shown). 

From these cumulative distributions and the KS-tests 
(Table [5]), the null hypothesis that the samples are drawn 
from the same population can clearly be rejected when com- 
paring the MM-only sample with those sources associated 
with both a methanol maser and radio continuum source 
(i.e. class MR). This suggests that the MM-only sample is 
not from the same parent distribution as sources with both 
a maser and radio continuum source for the mass parame- 
ter. Notably, the null hypothesis can also be rejected when 
comparing the MM-only sources with the combined star for- 
mation sample (class M -I- MR -I- R). This is not surprising 
considering that this sample is comprised of sources that 
have already been proven to be distinct from the MM-only 
sample (i.e. class MR). There is also a weak suggestion that 
the MM-only sample is not from the same distribution as 
sources with a radio continuum association (class R) for the 
mass. No distinctions could be discerned regarding the mass 
of the methanol maser sample (class M) . 

If we compare the cumulative mass distribution (Fig. [21 
middle) with the same plot produced in Paper I (Fig. 4) for 
an assumed temperature of 20 K, we find many similarities. 
Fig. 4 of Paper I showed that the MM-only sources were 
the least massive of the sample, followed by the methanol 
maser sources, the radio continuum sources, with the class 
MR sources the most massive in the sample. This result 
also holds true for the cumulative mass plot in this paper, 
which depicts a more robust estimate of the mass, taking 
into account the uncertainties on the temperature, rather 
than assuming a fixed one. 

The cumulative distribution of the combined star for- 
mation activity sources (classes M, MR and R) indicates 
that these sources are more massive on average than the 
MM-only sources. Median values for each sample are given 
in Table [4] which further corroborates this result. 

Figure [2] (bottom, left) shows the 2-D probability plot of 
the temperature and the mass. It is clear from this plot that 
these two parameters are highly correlated, which is hardly 
surprising given the nature of the modelling. The shading on 
the plot indicates the probability of occurrence for both the 
temperature and mass of the source in our sample, marginal- 
izing the contribution of the hot component for those sources 
where it is applicable. From this plot it is possible to ascer- 
tain the most likely combination of temperature and mass 
values for the sources in our sample. 

The 2-D probability plot corroborates the cumulative 



plot of temperature indicating the low probability of the 
sources in our sample having temperatures of <5K. The 
most likely temperature of the sources in our sample is be- 
tween 10 and 20 K, again confirming that our sources are 
cold, as well as the results of our earlier work in Paper I. 
The 2-D probability plot also indicates that the sources in 
our sample span a wide range of mass values, with the most 
probable mass roughly around 10* M© . The most likely com- 
bination of these parameters for our sources is a tempera- 
ture of 17 K and a mass of 5.7 X 10" Mq. This plot also il- 
lustrates that the very low and high temperatures (< 5 K 
and > 50 K respectively) correspond to very high and low 
masses (> 10^ M© and < 10 M0, respectively). Although, 
these combinations of parameters are plausible solutions 
to the SED analysis, they do not correspond to physically 
meaningful solutions, and the low, but non negligible, asso- 
ciated probabilities should be interpreted with care. 



4.4 Luminosity 

The cumulative distribution of the source luminosity is also 
presented in Fig. [2] (bottom, right). These plots, in addition 
to the KS tests, indicate that there is little difference, if 
any, between the different classes of sources in the sample 
in terms of their luminosity. Again, we attribute little sig- 
nificance to low luminosity values in the distribution, which 
are simply the tail-end values of the probabilities which ac- 
cordingly have a low probability of occurrence. Interestingly, 
the distributions for each of the different classes of source in 
the sample do not display the same shape or gradient. The 
mean and median values of the luminosity for each of the 
different classes in the sample are presented in Table 21 



5 DISCUSSION 

5.1 Massive star formation: luminosity, mass and 
SED modelling 

Identification and characterisation of young massive stars at 
all stages of their evolution, especially of the earliest stages, 
is essential in addressing and defining evolutionary scenarios 
for their formation. Among the meaningful physical param- 
eters, the mass and luminosity estimates of massive star- 
forming cores could play a pivotal role in their characterisa- 
tion and classification, possibly providing insight into their 
evolutionary status in a similar way that they do for low- 
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mass protostar classific ation from class to class I objects 
(e.g. lAndre et al.]|200Ci ). 

There are many diflterent data analysis techniques to de- 
termine the mass and luminosity of candidate young massive 
(proto)stars. Assuming a temperature, the mass estimate is 
usually derived fro m knowledge of the (sub) millimetre con - 
tinuum fluxes (e.g. iMotte fc Andrei I2OOII : iHill et ahllioosl) . 
Whilst assuming a temperature across a source sample is 
useful for first order approximations of the physical proper- 
ties (e.g. mass) of a so urce, in real ity all sources in a region 
iMotte fc Andrei I2OOII ) or survey (|Hill et al.ll2005l ) will not 
be at the same global temperature. We cautioned in Paper I 
that temperature assumption for mass determinations could 
lead to under- or over-estimations of the mass if the temper- 
ature was in fact cooler or warmer, respectively. It should be 
expected that different classes of source, such as methanol 
maser and radio continuum sources, which originate under 
different conditions, will also have different physical charac- 
teristics - including (but not limited to) their temperature 
and mass. For a large sample size, comprised of different 
classes of source, we therefore should not expect a global 
temperature to apply to a particular star formation com- 
plex or sample. 

SED modelling could, in contrast, provide us with useful 
estimates of the luminosity, temperature, and mass of star- 
for ming cores if the observational data are well constrained 
fe.g. lBurton et al.ll2004l : lMinier et al.ll2005h . In this instance, 
data are typically fitted with SED components that repro- 
duce the cold core emission in the far-infrared-millimetre do- 
main and the warm infrared emission that could be caused 
by a cluster of young stars. The shape of a SED, as well 
as the position of the SED majcima in the far-infrared- 
submillimetre and mid-infrared domains, can be indicative 
of evolutionary stage (e.g. iMinier et al.ii200f? ). Compiling a 
high quality well-constrain ed data set for sp ectral energy 
distribution modelling (cf. iMinier et al.l |2005| ) is, however, 
difficult, especially with the present lack of suitable high- 
resolution data covering the peak of the SED at ~ 100 fim. 
The extension of this method to large source numbers makes 
this an arduous task at best. This method is limited to small 
numbers of well-studied individual sources. 

Various radiative transfer models have also been devel- 
oped which build large grids of SED models for young stel- 
lar objects (YSOs) to which observational data can be com- 
pared a nd fitted. Howev e r, rad iative transfer models, such a s 
those of IWhitnev et all (|2003l ) and iRobitaiUe et all (|2006h . 
have been mainly developed for low mass (proto)stars. 
These authors are principally concerned with modelling mid- 
infrared emission from low-mass young stellar objects in rel- 
atively nearby star forming regions. Their conclusions rest 
on the characteristics of the warm, mid-infrared dust emis- 
sion surrounding these (proto)stars, and their orientations to 
our line of sight. Applying thes e SED fittin g techniques has 
recently been attempted bv.Mo linari et al. (2008). However, 
there are many caveats and limitations to SED modelling of 
young massive (proto) stars such as the presence of multiple 
sourc es and the confusion due to large distances (jRobitailld 
I2OO8I ). which need to be taken into consideration. 



5.2 Our approach: SED modelling with Bayesian 
inference method 

In our approach to SED analysis, we adopted a very simple 
description for the calculations of synthetic SEDs, rather 
than using detailed radiative transfer modelling. This ap- 
proach was chosen for two reasons: 

Detailed radiative transfer modelling requires a large 
number of free parameters (generally more than 10) which 
describe the spatial distribution of the dust grains and the 
source illumination, for example. When the available obser- 
vational data is limited, this kind of modelling results in 
degeneracies between each of these parameters, preventing 
firm constraints from being established. As we have only 
partial wavelength coverage of the spectral energy distribu- 
tion, only the characteristic parameters of our sources, such 
as the temperature, mass and luminosity, can be extracted 
in a quantitative way. We thus chose to restrict the number 
of free parameters for our SED modelling. Although provid- 
ing a less accurate description, this simple model provides 
a valid alternative to estimate the main parameters of the 
sources. 

In addition, only simple, analytical models allow sys- 
tematic sampling of the parameter space in a reasonable 
period of time (6.25 x 10® models in a four-dimensional pa- 
rameter space). This subsequently allows statistical stud- 
ies of a large number of models, via Bayesian inference, to 
be undertaken in order to determine the interplay between 
parameters and establish robust ranges of validity on the 
parameters. Such a work, especially when dealing with a 
large number of sources, is beyond current modelling capac- 
ities if radiative transfer models are used. For instance, a 
large numerica l effort of SED m odelling of YSOs has been 
undertaken by iRobitaille et al.l |2006 ) , whose current grid 
of models includes a total number of 200 000 SED models, 
sampling a 14-dimensional parameter space. This limited 
number of models cannot provide a systematic sampling of 
the parameter space and the authors needed to bias the 
explored parameter space to accommodate the computing 
time required. As a result, trends and correlations between 
parameters are difficult, if not impossi ble, to apprehend with 
such a grid of models (see iRobitaillel (j2008) for a more de- 
tailed description of the limitations in the use of their grid 
of models to estimate parameters). 

5.3 Interpreting the distinction between mass and 
luminosity of our source samples 

Using the first application of the Bayesian inference method 
of SED modelling to massive star formation, we have es- 
timated a range of values for the temperature, mass and 
luminosity for sources from the SIMBA survey of iHiU et al.l 
(2005). We have performed a statistically robust analysis of 
each of these parameters across the different classes of source 
in the sample, with the aim of characterising and classifying 
them. 

If we attribute some significance to the small distinction 
between the cumulative distributions of the various object 
classes, MM-only cores to UC Hii regions (Fig. [2] and Table 
O, two conclusions can be reached. Firstly, the cumulative 
distribution of the methanol maser sources, for each param- 
eter, is very much representative of the entire sample. This 
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Figure 3. Mass-Luminosity diagram for the different classes of 
source in the sample. Each of the points represents the median 
value (i.e. the values at which the cumulative probability distri- 
bution is 0.5) of the mass and luminosity of the individual sources 
(green cross: mm, red plus sign: m, blue triangles:mr and pink cir- 
cles: r). The contour levels represent the regions enclosing 68 per 
cent of the Bayesian probability for a given class of sources (green 
full line: mm, red dashed line: m, blue dot-dash line: mr and pink 
dot-dot-dash line: r). The full and dashed straight lines depict the 
M = L and M = L"-^ relations, respectively. 



result (Fig. [2]) coupled with the fact that methanol masers 
are exclusive signatures of high mass star formation con- 
firms that our sample is representative of high mass star 
formation. In addition, the MM-only cores are less massive 
on average than cores with a methanol maser and/or radio 
continuum association and hence those sources known to 
support massive star formation. This is in agreement with 
Paper I, despite the different approaches to mass determina- 
tion. The cumulative distribution plots of the temperature 
and luminosity (Fig. [2J display only small distinctions be- 
tween the different classes of source in the sample, which 
is more the case for the luminosity than the temperature. 
Notably, there is no significant difference between the MM- 
only sample and the combined star formation sources, i.e 
those sources with a methanol maser and/or radio contin- 
uum source (class M, MR and R) . This result suggests that 
the MM-only sample has similar characteristics as sources 
with known star formation sites. However, the KS-tests in- 
dicate that the MM-only cores are not from the same parent 
population as sources with both a methanol maser and ra- 
dio continuum association (class MR) with respect to their 
masses. 

Figure [3] introduces a luminosity vs. mass diagram 
(hereafter, M — L diagram) for our sample of cores. A similar 
diagram was proposed for cla ss and class I low-mass pro- 
tostars by I Andre et al.l l|2000l ). who presented the protostel- 
lar envelope mass vs. the bolometric luminosity for individ- 
ual protostars or multiple protostar systems. Comparatively, 
in this work, the derived masses in Table [5] overestimate 
any protostellar envelope mass because the (sub) millimetre 
emission fluxes are integrated over more than a FWHM 
beam. The masses reported in Table |5] almost certainly in- 
clude contributions from both the core and from diffuse ex- 
tended gas, whilst the luminosities are rough estimates of 



the source bolometric luminosity. Consequently, if a domi- 
nant source is responsible for both the luminosity and the 
mass of the core, then the M — L diagram could prove a 
useful diagnostic of the evolutionary status of this dominant 
object. 

Figure plots the median value (the value correspond- 
ing to 0.5 on the cumulative distributions) of mass and 
luminosity for each of the sources in the sample, as well 
as the 68 per cent Bayesian probability of occurrence for 
each source class (the latter is explored in the following sec- 
tion). The masses and luminosities span 10 - 10** M0 and 
10^ — 1O®L0, respectively. These ranges a r e com parable to 
those presented in Fig. 9 of iMolinari et all (|2008|). although 
our ob jects have higher masses on average. Molinari et al.l 
(|2008t ) ' undertook a method of approach different to our 
own, in which they assigned millimetre emission fluxes to 
the massive star-forming objects by removing emission from 
the diffuse clumpy medium. Their resultant masses are con- 
sequently lower for their sources compared with this work. 
They also fitted their SEP to a grid of mo dels that were 
computed according to I Whitney et af] (|2003l ). 

In Fig. [3] the MM-only cores appear to segregate into 
two groups with the division loosely defined by the line of 
M = L'''^. The first group of MM-only sources lies in the 
same region of the plot as the sources with known star for- 
mation activity (methanol maser and/or radio continuum 
sources) i.e. classes M, MR and R. The other group of MM- 
only sources exhibit lower masses or higher luminosities and 
are clearly separated from class M, MR and R sources. 

In summary, the cumulative distribution plots produced 
from our SED models as well as the KS-tests and the M — L 
diagram reveal that some MM-only sources are distinct from 
those sources with a radio continuum association in terms of 
the mass parameter. These results are consistent with Paper 
I. In Paper I, it was also revealed that the MM-only cores 
have the smallest radii of the sample and at least 45 per 
cent are without mid-infrared MSX emission. The charac- 
teristics of some MM-only cores — the least luminous and in- 
frared quiet ones — are consistent with less evolved examples 
of massive star formation. However, we cannot exclude the 
possibility that some MM-only cores may be dominated by 
luminous intermediate-mass protostars, hence less massive 
but still relatively luminous cores. Finally, some fraction of 
the MM-only cores may even represent t he 'fa iled' cores hy- 
pothesised by 3^azquez-Scmadeni ct al] (|2005l ). More likely, 
the MM-only cores represent various classes of prestellar and 
protostellar sources. 



5.4 Significance of luminosity vs. mass diagram 

Alternatively, if we do not attribute any significance to 
the small distinction between the various source samples in 
Fig. [2] these results may simply reflect the limited informa- 
tion that can be extracted from our SED analysis, especially 
considering the poorly constrained data in the far-infrared 
regime. There are many limitations to the interpretations 
of spectral energy distributions and M — L diagrams for 
high-mass star forming cores. For instance, the luminosity 
estimates are poorly constrained in the far-infrared domain 
and the angular resolution of SEST at 1.2 mm is insufficient 
to excl ude the possibili t y of multiple source components - 
indeed [Longmore et al.l 1 200d ) have detected a few massive 
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infrared protostars or young stellar objects within 6,000 AU, 
corresponding to 24 arcsec at 4kpc (i.e. with our current 
data, we would not expect to see this). 

To account for this. Figure [3] also presents the 2-D 
Bayesian probability contours enclosing the 68 per cent 
Bayesian probability of occurrence for both the mass and 
the luminosity. It can be seen that the mass and the lumi- 
nosity encompass a much larger range of values than de- 
picted by the median values on the plot (see section [5.3(1 . 
The MM-only cores could occupy the range 100 Mq and 
10*^ Lq or 10'' Mo and 10 Lq of Fig. [3] with equal probabil- 
ity. That is, these results suggest that the MM-only cores 
could either be star-forming regions in quite an advanced 
stage hosting embedded, luminous stars (10^ L©) or alter- 
natively very massive quiescent clouds (10 ). Moreover, 
we observe a similarly wide range of probable luminosity 
and mass for the UC Hii regions and the methanol maser 
sources in our sample, despite the fact that these sources 
are well identified at high angular resolution. The Bayesian 
inference method demonstrates that SED fitting, with the 
currently available far-infrared data sets, cannot provide us 
with reliable evolutionary tracks in the M — L diagram for 
high-mass star formation. 



6 CONCLUSIONS 

We have performed spectral energy distribution analysis 
for 227 of th e 405 sources detected in the SIMBA survey 
of iHill et al.l (^2005). Using the Bayesian inference method 
of analysis we have determined a range of suitable values, 
with associated probabilities, for each of the parameters of 
temperature, mass and luminosity. Each of these parame- 
ters have been analysed with respect to the different type 
of source in the sample and hence their associations (or 
lack thereof) with methanol maser and/or radio continuum 
sources. The cumulative distribution plot of the mass for the 
different source classes is consistent with our earlier work. 

If we attribute little significance to the class distinction 
for each parameter of mass, luminosity and temperature, 
then the MM-only cores have the same characteristics as 
sources with known star formation activity (classes M, MR 
and R), yet they display no overt signs of star formation. 
Following this, the MM-only cores are excellent candidates 
for early stage protostars or massive young stellar objects. 

Attributing significance to the marginal distinctions be- 
tween the MM-only sources and those sources with methanol 
maser and/or radio continuum associations for each of the 
temperature and luminosity, and factoring in the results 
of the mass, radius and lack of mid-infrared associations 
for almost half, then we can interpret that the MM-only 
core is a younger, smaller and less evolved example of mas- 
sive star formation. That is, they represent an evolutionary 
stage of massive star formation, prior to the development of 
methanol maser emission and are thus indicative of the earli- 
est stages of massive star evolution. However, we cannot ex- 
clude the possibility that the MM-only cores are examples of 
'failed' cores or instead will support intermediate mass star 
formation. Alternatively, the MM-only core could comprise 
a cross-section of sources supporting both arguments. Spec- 
tral line observations (e.g. of turbulent linewidths and/or 
chemical state) of these MM-only cores are necessary in or- 



der to determine which of them, if any are forming massive 
stars. 

It is clear from this work, that SED modelling is heavily 
reliant upon well-constrained and robust data which is well- 
sampled in wavelength space. There is a clear dependency of 
the stringency of the fit upon the quality of the data. From 
our Bayesian inference SED analysis, it is clear that in the 
absence of reliable far-infrared data, which would serve to 
constrain the peak of the SED, it is not possible to draw 
reliable evolutionary tracks in the mass versus luminosity 
diagram of high-mass star formation. Future observations 
with the Herschel Space Observatory will provide greater 
constraints in the crucial far-infrared/submillimetre regimes 
for SED modelling. 
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